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HYPERVELOCITY BINARY STARS: SMOKING GUN OF MASSIVE BINARY BLACK HOLES 

YOUJUN Lu\ QlNGJUAN Yu\ & D. N. C. LlN^'^ 
Draft version February 1, 2008 

ABSTRACT 

The hypervelocity stars recently found in the Galactic halo are expelled from the Galactic center 
through interactions between binary stars and the central massive black hole or between single stars 
and a hypothetical massive binary black hole. In this paper, we demonstrate that binary stars can be 
ejected out of the Galactic center with velocities up to 10'^ km s~^, while preserving their integrity, 
through interactions with a massive binary black hole. Binary stars are unlikely to attain such high 
velocities via scattering by a single massive black hole or through any other mechanisms. Based on 
the above theoretical prediction, we propose a search for binary systems among the hypervelocity 
stars. Discovery of hypervelocity binary stars, even one, is a definitive evidence of the existence of a 
massive binary black hole in the Galactic center. 

Subject headings: black holes physics - Galaxyicenter - stellar dynamics 



1. INTRODUCTION 

Evidence of a massive black hole (MBH) is securely 
estab l ished in the Galac tic center (GC) ( Schodcl et aL, 
20021: iGhez et al.l l2003l [20051 : lEisenhauer et all I2005h . 
Recent discoveries of young sta.rs in the sub-parsec re- 
gion around the GC (jSandersI [Tool IGhez et al.l [200l 
and h ypervelocity stars (HVSs: first recognized by 'Hiiy 
19881) in the Ga l actic halo (Brown ct al. 2005, 2006a^b; 



Edelmann et all 120051 : iHirsch et all 120051) have stimu- 



lated the hypothesis of a second MBH or an intermediate- 
mass BH (IMBH) which may be common in galactic nu- 
clei if they are assembled through mergers of smaller 
galaxies and/or star clusters with central MBHs or 
IMBHs (Bcgcl man et"all [l980l: I^IMl IVolonteri et all 
[2003 ). IMBHs can transport young stars to their cur- 
rent locations (jHansen fc Milosavlievid [2003f ) where in- 
situ star formation is suppressed by the strong tidal force 
of the MBH. Through close encounters, a massive bi- 
nary black hol e (BBH) may also eje c t nearby stars to 



becoine HVSs (IYu fc Tremaine 



[20051: [Le^[20M iSesana et al 



2003t iBaumgardt et al.l 
2006D . Direct observa- 



tional evidences for a secondary MBH are difficult to es- 
tablish because the semimajor axis of the BBH could be 
well within the orbit of the most central 'S' star s which 
have longer- 1 han- a- dec ade orbital periods ((Schodel et al.l 
[200l[caleri~aI][2?Ki3h . 

With a simple scaling argument in § [2] and a series of 
numerical calculations in § [3| and § |31 we demonstrate 
that hypervelocity binary stars (HVBSs) with velocities 
up to 10'^ km s~^ can be ejected out only by their dy- 
namical interactions with a massive BBH. We propose 
a search for HVBSs among HVSs. Discovery of HVBSs, 
even one, in the Galactic halo will be a definitive evidence 
for the existence of a massive BBH in the GC. 

2. DOMAINS OF RELEVANT DYNAMICAL PROCESSES 

We characterize the dynamics of stars around a MBH 
in galactic centers with the following critical radii: (1) 
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the Schwarzschild radius is rsch = 2GMbh/c^ = 3.4 x 
10~''mBH pc, where tobh is the MBH's mass Mbh 
normalized b y that of the MBH in the GC (in unit o f 
3.6 X IO^Mq] IGhez et al.l [20051 : lEisenhauer et all 120051 ). 
G is the gravitational constant and c is the speed of 
light; (2) the tidal radius for a single star with mass 

1/3 

and radius is r^;^ — {rf'M-Qii/'m^^ = 3.5 x 
lQ-^{TfmYi]i)'^/'^{MQ/m^)'^l'^{r^/ Rq) pc, where a con- 
stant 77 = 2.21 and 0.844 for a homogeneous, incom- 
ressible body and an n = 3 polytrope res pectively 
Sridhar fc Tremaind [l99^ iDiener et al.lll995[ ): (3) the 
tidal breakup radius for a binary star (with semima- 
jor axis ttb and component masses mi and TO2, re- 
spectively) is Tj^id = ab[2MBH/(n^i + ^2)]^/^ = 7.4 x 

10"''^"^BH(ab/0.1 AU)[2M0/(mi ma)]!/^ pc; (4) the 
radius of the sphere of influence of the MBH is OeS = 
GA/bh/ct'^ ~ 2.3 pc, where a is the ID velocity dis- 
persion of the spheroidal component of the galaxy and 
the Mbh — relation (Mbh cx: cr^°^; I Tremaine et al.l 
l2002f ) is used to replace a by Mbh; and (5) a mas- 
sive BBH becomes "hard" when its semimajor axis 
Qbbh ^ o,h, where the critical value Oh = G'AfBH.2/4cr^ ~ 
5.9 x 10-3 mO|(j./0.01) pc, ly = Mbh, 2 /Mbh, Mbh = 
Mbh.i + Mbh,2 (Mbh,i > Mbh, 2), and Mbh,i and 
Mbh, 2 are the mass es of the primary and secondary 
MBHs, respectively (jOuinlanl [l996) . BBHs, if exist 
in nearby normal galactic centers, would spend most 
fraction of their lifespan a t agBH O.Oloh — lah if 
lO^M© < Mbh < lO^M© (fYul[2003 l. so that for com- 
panions around these MBHs with > 10"^, we generally 
have rsch < r^-^ < r^id ^ "bbh < Oh < Ocff (see Fig. [T|). 

With these definitions, the vicinity of a MBH can 
be partitioned by various dominant physical processes. 
Most of the stars within r < Ocff have negative en- 
ergy and are bound to the central MBH and their kine- 
matic properties have been u sed to determine the MBH 
mass (jRichstone et al.l 119981 ). Stars that enter into the 
region r < rsch (if ''tid < ^sch) are swallowed as a 
whole by the MBH. Stars that enter into the region with 
''Sch < f ^ ^tid tidally disrupted. A fraction of the 
disrupted star is accreted by the MBH through a disk- 
like flow with outbursts of enormous radiation and the 
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other fraction m ay be ejected out of galactic nuclei (|Hillsl 
[19751: lReeslfT988h . 

Here, we are primarily interested in the origin of 
HVSs. A star initially unbound or weakly bound to a 
MBH will attain a velocity v ~ ^ GAlBn/r ~ 5.5 x 
10"^ km s^^mg'^j^(l mpc/r)^/^ upon moving within a dis- 
tance r from the MBH, where mpc = 10~'^ pc. If the star 
then ensures a velocity perturbation dv <^ v with an in- 
crease in the specific energy of the star v6v which is much 
larger than the magnitude of its initial specific energy 
(|w2/2- GMbh/^I), it would escape from the MBH with 
velocity ^ \/2v5v. For example, a binary star that enters 
into the region r^j^j < r < rj'jj is likely to be tidally broken 
up, during which episode, each star would receive a 5v 
comparable to its binary orbital velocity. Consequently, 
one component of the binary star would be ejected out 
of the system with hyper velocity up to 10"^ km s~^ and 
the other w ould become more tightly bound to the BH 
(|Hillslll988l ). Around a hard BBH with obbh ^ flh, most 
low-angular-momcntum single stars that enter into the 
region r^j^ < r < aBBH can also be ejected by the BBH. 
The r.m.s. of the velocities of the ejected stars at infinity 
is given by 

vTj ^ y^2i^GAfBH.a^BH,2/(MBHaBBH) 

^900 km s-im5^-25(l _ i.)i/2(o.iah/aBBH)'/', (1) 

(see eq. 17 in rYu|[2002l ). where K ~ 1.6 is a constant 
and the Mbh ~ relation obtained by Tremaine et at 
(2002) is adopted. The main portion of Sv is due to 
the impulse induced by the secondary MBH F5t, where 
F ~ GMBH,2/aBBH the force per unit mass from the 
secondary BH and St ~ GMbh/abbh interaction 
time. Provided that obbh is substantially larger than 
rj'jjj, there is a large probability for a binary star's peri- 
apse not to enter into the region r < r^^^ (with respect 
to either component of the BBHs) before it is ejected. 
Consequently, binary stars which enter into the region 
r^^^ < r < cbbh niay be ejected, with large velocities 
while preserving their systems' integrity, as HVBSs. 

3. NUMERICAL METHOD 

In order to determine the probability of HVBS for- 
mation, we simulate the interaction of binary stars with 
BBHs. Instead of carrying out a comprehensive series of 
prohibitively complex four-body problem, we note that 
the motion of a binary star with masses mi and m2 in 
the potential of a BBH may be greatly simplified in the 
limit of mi, 7712 ^ A/bh.i, AfBH,2 because the motion of 
the BBH is not affected by the binary star and can be 
described in terms of a two-body problem. We set the 
center of mass of the BBH is set to be at rest at the origin 
of the coordinate system. Each component of the binary 
star moves in the potential of the rotating BBH and that 
of its stellar companion (i.e., Gmj/lfj—ril—Gmjri-rj/r^, 
i ^ j, i,j = 1, 2, where rt and rj are the position vectors 
of the two components). Here Newtonian mechanics is 
applied. The initial velocity of the center of mass of the 
binary star at infinity is set to be the velocity dispersion 
a of the galaxy. We choose its impact direction randomly. 
The impact parameter b is chosen so that the pericenter 
distance rp of its corresponding Keplerian orbit [which is 



related with b through b = rp(l + 2GMBii/ {cr^rp)y^^] is 
in the range from CIobbh (or COIobbh) to 2aBBH- The 
angular momentum vector and the initial phase of the 
binary star's orbits about their center of mass are chosen 
randomly. For each set of initial conditions, we inte- 
grate a system of differential equations with an explicit 
5(4) order Runge-Kutta scheme. The center of mass of 
the binary star is analytically extrapolated from r — oo 
to IOOcbbh along a Keplerian orbit about a point mass 
A^BH where the numerical integration starts. 

Some integrations are time consuming because the bi- 
nary is captured into weakly-bound orbits to the BBH 
and makes many revolutions before it is expelled. We 
adopt similar approximations to expedite t he nu merical 
integrations as that introduced bv lQuinlaiil {l99^ for the 
studies of interactions between a single star and a BBH. 

The fate of the interaction of a binary star with a mas- 
sive BBH mainly falls into two cases: (1) the binary star 
is ejected as a whole from the BBH with some changes in 
its intrinsic orbital semimajor axes and eccentricities; or 
(2) the binary star is dissociated by its interactions with 
the BBH and both dissociated components are ejected 
with hyper velocities. 

4. NUMERICAL RESULTS AND DISCUSSIONS 

We present numerical results of the interaction of a 
binary star with a massive BBH in Figures [2] and [3l 
We assume various values for the semimajor axes of bi- 
nary stars (with Cb — 0.1 AU and 0.3 AU; note that 
here 0.3 AU is taken as an upper limit because bina- 
ries with larger ab are likely to be disrupted by encoun- 
ters with single stars in the GC within a Hubble time, 
lYu fc Tremai^l2003 ^ and for the BBH semimajor axes: 
OBBH = O.loh = 0.37,5.88,0.48 mpc [= (4.9,23,7.2) r"^-^ 
for the Milky Way, M31, and M32, respectively] and 
OBBH = arj'ij = 0.15 mpc (for the Milky Way). The 
component masses of the binary star are set to be 777 1 = 
77i2 — IMq and our results are not modified by other 
values of binary masses. The BBH mass ratio is set 
to he ly — 0.01 unless otherwise stated. For simplicity, 
the BBH is set to be on a circular orbit and the initial 
intrinsic orbit of binary stars are also circular. Figure 
[^Ja) shows that the probability of the binary star being 
ejected away as a whole from the BBH is generally larger 
than 50% if rp > r^j^j and that binary stars are more likely 
to be dissociated with decreasing pericenter distances rp. 
Given the same BBH mass and semimajor axis, binary 
stars with larger Ob's are more likely to be dissociated 
(see the thick solid line and dot-dashed line). If we keep 
obbh = O.lflh but set — 0.001 for the assumed BBH 
in the GC, the probability of ejecting HVBSs is smaller 
(down to 20% at rp ~ obbh) since cbbh is smaller than 
the tidal radius of the primary MBH r^^^ and most of the 
binary stars are tidally broken up. 

Figure [3] shows the dependence of the ejection veloc- 
ities of HVBSs on the pericenter distance rp, which is 
generally consistent with the dependence of the ejection 
velocities of single stars through three-body interactions 
with the BBH (thin solid line). The can attain val- 
ues up to 10^ km s~^ at rp ^ obbh = O.lah and it re- 
mains roughly constant with decreasing rp but decreases 
with increasing rp. The ejection velocities of HVBSs de- 
pend on the BBH mass and semimajor axis (see eq. [1} . 
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Fig. 1. — Various characteristic radii in galactic centers with 
massive black holes: (1) the Schwarzschild radius (dashed line); (2) 
the tidal radius for a single star with solar mass and solar radius 
(dotted line); (3) the tidal breakup radius for a binary star with 
semimajor axis a^, = 0.1 AU (long-dashed line; r] = 2.21 is adopted 
here, see text); (4) the radius of the sphere of influence of central 
MBH (dot-dashed line); and (5) the transition radius for a hard 
binary black hole with mass ratio of i/ = 0.01 (thick solid line), 
the thin solid line represents aBBH = O-lajj. The characteristic 
radii for the MBHs in the Milky Way (with Mbh = 3.6 X IO'^Mq 
and a = 103 km s"!), M31 (1.4 X W^Mq, 160 km s"!), and 
M32 (2.5 X IO'^Mq, 75 km s"!) are labeled as circles, squares and 
triangles, respectively. 




Fig. 2. — Probability for ejecting hypervelocity binary stars, 
(a) Fraction of undisrupted binary stars that are expelled after 
their interactions with a hypothetical BBH. The BBH mass ra- 
tio is set to V = 0.01. The thick lines are for the BH mass in 
the Milky Way with different model parameters for the BBH and 
binary star's semimajor axes. The thin long-dashed line is for 
the BH mass in M31 and the thin dotted line is for M32. The 
circles/triangles/squares mark the radii rj^^ for each cases, (b) 
Fraction of binary stars that are expelled away as a whole with 
hypervelocity t;^ > 900 km s~^ after their interactions with the 
assumed BBH in galactic centers. The line types and the solid 
points have the same meanings as those in (a). As a reference, 
the thin solid line represents the fraction of single stars that are 
expelled with v > 900 km s~^ after their three-body interactions 
with an identical BBH as that represented by the thick solid line. 

With all else being equal, the ejection velocities increases 
with increasing BH mass (see the thin lines for M31 and 
M32, respectively, in Fig. [31); and decrease with increas- 
ing BBH semimajor axes (see thick short-dashed line and 
thick solid line). 

Figure [2Ib) shows the probability of the binary star 
being ejected away as a whole with hypervelocity > 
900 km s^^ from the BBH has a peak value of 10%- 
50% at Tp ~ (0.1-1.5)aBBH and it decreases for both 
larger and smaller values of Vp. (Due to the decelera- 
tion in the Galactic potential, an ejected binary star, 
HVBS, with vf^ = 900 km s'^ would attain an asymp- 
totic velocity 700 km s^^ at the Galactocentric distance 
50 kpc, which is corresponding to the observed radial 
velocities of HVSs in the Galactic halo.) The decline of 
ejection probability of binary stars with high velocities at 
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Fig. 3. — The r.m.s. velocity of the ejected binary stars at infinity 
versus the pericenter distances of their initial injection orbits. The 
line types have the same meaning as those in Fig. |2] 

large rp is consistent with that for single stars (see thin 
solid line). But the rarity of ejecting HVBSs at small 
Tp is mainly due to the dissociation of binary stars by 
the Hills mechanism as indicated in Figure [SJa). (Note 
that the encounters with a BBH at small Vp require low- 
angular-momentum orbits. After the depletion of such 
an initial population, the likelihood to eject HVBSs at 
small rp declines because, through two-body relaxation, 
the remaining binary stars would slowly lose their angu- 
lar momenta to other stars and reach the vicinity of the 
BBH preferentially with large rp.) Our calculations show 
that for a larger BBH semimajor axis obbh = 0.3ah in 
the GC, the peak probability is about 5%-10%, which 
is still non-negligible. The ratio of number densities be- 
tween the HVBSs and the single HVSs contains informa- 
tion on the fraction of binary stars in the GC, the BBH 
mass ratio and semimajor axis. If a BBH with u 0.01 
and flBBH ~ O.lflh exists in the GC, according to Figure 
[2ljb), the probability of binary stars being ejected with 
> 900 km s^i would be at least about 1/2 - 1/3 of 
the probability of single stars being ejected with such hy- 
pervelocities provided rp > r^^^. If the density of binary 
stars with Ob ~ (0.1 — 0.3) AU is about 10% of that of 
the single stars in the GC, the fraction of HVBSs among 
the HVSs would then be at least about 3% — 5%. 

The ejection rate of HVSs by an assumed BBH (with 
f = 0.01 and obbh = 0.5 mpc) in the GC would be 
~ 10~^ yr~^ if the low-angular-momentum stars initially 
in the loss cone are depleted and the loss cone is re- 
filled by two-body relaxation processes (|Yu fc Tremaind 
l2003f ). If the fraction of binary stars with ab ^ 0.3 AU 
in the GC is about 10% of total stars, the ejection rate 
of HVBSs would be > 5 X 10~^ yr~^. By apply i ng th e 
same method that is adopted bv lYu fc Tremaing (|2003f) . 
we obtain ejection rates for M32 (with i/ = 0.01 and 
obbh = 0.48 mpc), which is similar to that for the 
Milky Way But for M31 (with ly = 0.01 and obbh = 

^ The hardening timescale of the BBH due to three-body inter- 
actions with stars passing through its vicinity can be > 10^ yr, and 
the gravitational radiation is the dominant mechanism to make the 
B BH lose energy with a timescale of 10* yr. According to Figure 2 
in lYu fc Tremaind ll200a ). the BBH should have f > 10"^ in order 
to eject stars with ti^ > 700 km (which corresponds to a ve- 
locity > 400 km s~^ at the Galactocentric distance 50 kpc) and to 
have the present accretion rate of BH companions is not unusually 
large compared to the average over the lifetime of the Galaxy. 



4 



Lu, Yu, & Lin 



5.9 mpc), we obtain lower ejection rates (~ 10 ^ yr ^ 
for HVSs and ~ lO-^-lQ-^ yr"! for HVBSs). For M31, 
we adopt the surface densi ty distribution estimated by 
iKormendv fc Bender! ()1999l . eq. 1) and we neglect the de- 
parture of ste llar distribution from spherical symmetry 
([Bender et alJ 12005 ) as for the other cases. Asphericity, 
if significant, would enhance the ejection rates of HVSs 
and HVBSs. 

Several HVSs have recen tly been discovered throu gh 
spectroscopic observations (Brown ct al.' l2005l l2006al fbl: 
lEdelmann et al. 2005; Hirs.ch_ct_al,, ,200^. The pres- 
ence of any HVBSs among this sample is difficult 
to be spatially resolved. But the spectrum of the 
more conspicuous primary would exhibit measurable (~ 
100 km s~^sini ) periodic variations on the timescale 
of a few to several ten days. Transits may also intro- 
duce observable periodic light curve modulations. Non- 
detection of HVBSs does not rule out the existence of 
a BBH in the GC since their absence can be attributed 
to: (1) a low rate of binary-stars injection into the re- 
gion rji'jj < r < obbh; (2) a relative compact BBH with 
obbh < f'Jid; or (3) a relatively soft BBH which is unable 
to eject stars with sufficiently high velocities. 

In contrast, the detection of any HVBSs unequivocally 
proves the existence of BBH because the production of 
HVBSs by any other mechanism appears to be challeng- 
ing and unlikely. (1) The Hills mechanism cannot lead to 
the production of HVBSs because it requires the breakup 
of a binary star by a single MBH. (2) In principle, tidal 
disruption of a hierarchical triple star (a close binary plus 
a less tightly bound tertiary companion) by the central 
MBH may lead to the capture of the tertiary by the MBH 
and the ejection of the close binary as a HVBS. However, 
the rate of ejecting HVBSs by this mechanism is negli- 
gible if the following factors are taken into account, (a) 
The fraction of the binaries with additional companions 
(triple, quadruple, or quint uple) is about 60% of t he bi- 
naries observed locally bv iTokovinin et aT] (|2006f l. (b) 
All known tertiary companions around the binaries have 
periods longer than 2 yr (with semimajor axis > 2 AU) 
and >90% of the tertiary companion have orbital peri- 
ods ranging from 5 yr to 10^ yr. (c) The typical ejec- 
tion speed of the close binary due to the breakup of 
a hierarchical triple with semimajor axis 2 AU is only 
~ 400 km s~^. If the ejection speed follows a Gaussian 
distribution with a dispersion which is about 20% of the 
average speed (the dispersion is mainly due to the dif- 
ferent orbital phases and orientations of t he binary orbit 
relati ve to the MBH; e.g., see Fig. 1 in iBromlev et al.] 
I2006D . stars with velocity higher than 700 km s~^ or 

900 km s~^ would be beyond 3a or 5tT from the mean, 
which corresponds to small probabilities of 1.3 x 10~'^ or 
3 X 10~^. Thus, the ejection rate of HVBSs due to the 
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tidal breakup of hierarchical triple stars is only a fac- 
tor of 0.6 X 0.1 X 1.3 X IQ-^ - 8 X 10"^ or 2 x 10"** of 
the number of HVSs due to Hills mechanism. We note 
that the argument above has some uncertainties because 
the estimate derived for the solar neighborhood is uncer- 
tain and it is not clear whether the triple distribution 
at the Galactic center is the same as that in the solar 
neighborhood. (3) HVSs may also be ejected as stars are 
scattered off a cluster of stellar-mass black holes (SBHs) 
orbiting a central MBH. In principle, the ejected stars 
can be binaries. But, the perturbation on the velocity 
of ejection star(s), Sv is primarily due to close encoun- 
ters with slightly - more-massive SBHs at a distance about 
the solar radius (jO'Learv fc Loebll2006D . Such close en- 
counters can lead to either the disruption of the binary 
system or an exchange in which the SBH captures one 
component of the initial binary while the other is ejected 
as a single star. The newly formed star-SBH binary can- 
not be ejected out of the GC as a HVBS since the SBHs 
are closely bound to the central MBH with a substan- 
tially negative specific energy {\E\ 3> v6v because usu- 
ally Sv <C v). Whereas the detection of a H VS provides 
support for the exi stence of a MBH in the GC (lHillslll988l : 
iBrown et aLll2005f) . the detection of a HVBS will be the 
smoking gun for the existence a massive BBH in the GC 
(or M31, M32). 

Our numerical calculations show that both Ob and the 
eccentricities of the binary may be strongly modified dur- 
ing their interactions with the BBH. The eccentricity of 
the binary-stars' orbits can be excited up to 0.9-1, which 
may lead to the merger of the two components of the 
binary star with small relative velocity. Therefore, the 
interactions between binary stars on bound orbits and 
a BBH may provide a channel for the formation of the 
S-stars discovered in the Galactic center. It is plausible 
that both the formation of S-stars and the ejection of 
HVSs are the by-products during the process of a MBH 
inspiralling inward towards the primary MBH. 

The velocity of some binary stars ejected by interac- 
tions with a BBH can be smaller than the escape veloc- 
ity of the Galactic halo and these stars are bound to the 
Galaxy. Such a binary may already be contained in ex- 
isting observational data. For example, Scorpius X-1, a 
X-ray binary, is found to have a velocity of ~ 480 km s~^ 
with its past pe rigalactic distance of ~ 500 pc to the GC 
(jMirabel fc Rod rigucs 2003). This object could be ex- 
plained as a binary formed in the GC and ejected out by 
a BBH from there, though other explanations, such as a 
natal kick of supernovae explosion, are possible. 
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